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ABSTRACT 

ANALYSIS  OP  INELASTIC  ALPHA  PARTICLE  SCATTERING 
USING  DIRECT  INTERACTION  THEORY 

by 


Herbert  Prank  Swanson  Jr. 

The  theory  of  the  direct  interaction  is  discussed 
and  its  applicability  to  inelastic  alpha  particle 
scattering  to  low  lying  states  is  shown.   Two  results 
of  this  theory  are  developed;   (1)  the  diffraction  model 
of  Blair  and  (2)  the  distorted  wave  Born  approximation, 
where  the  phenomeno logical  approximation  to  the  optical 
wave  function  suggested., by  McCarthy  and  Pursey  is  used. 
Fits  to  Ar3b(  o4,  <*')Ar^b""'  data  are  shown  for  32.8  and 
41.0  Mev(Lab)  alpha  particles.   In  the  diffraction  model 
analysis,  RQ  =»  6.45  Fermi s  for  both  energies. 
/3>RQ  =  .82  Fermi s  for  the  lower  energy  data  and  .71 
Fermis  for  the  higher  energy  data.  /4l?Ro  =*  .67  and  .56 
Fermis  for  the  lower  and  higher  energy,  respectively. 
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CHAPTER  I 
THE  DIRECT  INTERACTION 

There  are  two  basic  mechanisms  by  which  a  nuclear 
reaction  can  proceed;  (1)  the  direct  interaction,  and 
(2)  the  formation  of  a  compound  nucleus.  The  theoretical 
angular  distribution  predicted  by-  each  mechanism  is 
different.  To  see  the  nature  of  these  two  processes, 
consider  the  incident  particle  to  be  a  nucleon.  Thus 
the  independent-particle  shell  model  is  a  good 
approximation  to  the  nuclear  structure.   Here  the  many 
body  problem  is  replaced  by  an  effective  single  particle 
potential  which  describes  the  force  acting  on  each 
particle,  this  potential  being  different  for  different 
shell  model  configurations. 

As  a  particle  enters  the  nucleus,  it  experiences 
a  potential  described  by  the  other  A  nucleons  in  their 
ground  state. configuration.  The  projectile  may  then 
simply  enter  the  nucleus,  be  deflected  by  the  nuclear 
potential,  and  emerge  again  at  a  different  angle  but 
with  the  same  energy  in  the  center  of  .mass  system. 
This  is  called  direct  elastic  scattering.   If  a 
reaction  ensues,  the  energy  transferred  to  the  nucleus 
is  determined  by  the  interaction  of  the  incoming 
nucleon  with  one  in  some  shell  configuration  and  by 


the  nuclear  density.   This  later  fact  modifies  the 
problem  as  neither  of  the  colliding  particles  can  have 
an  energy  below  that  of  the  ground  state  energy  of  the 
nucleus.   This  is  because  the  hole  left  by  the  struck 
nucleon  is  the  only  level  available  below  the  ground 
state  energy.   The  case  where  the  incident  particle 
replaces  the  struck  nucleon  is  called  exchange  scattering 
and  will  not  be  considered  here.   Those  nucleons  excited 
to  high  level  states  should  be  described  by  a  different 
potential  than  those  excited  to  low  lying  states, 
hence  the  effective  potential  is  a  function  of  the 
bombarding  energy.   If  either  of  the  colliding  nucleons 
has  an  energy  greater  than  its  separation  energy,  it 
may  leave  the  nucleus  without  further  reaction  save 
deflection  by  the  nuclear  potential.   This  is  described 
as  a  direct  interaction  and  usually  occurs  in  a  time 
comparable  with  the  time  it  takes  the  projectile  to 
traverse  the  nucleus.   If,  on  the  other  hand,  neither 
of  the  nuclei  has  energy  sufficient  for  separation, 
both  will  undergo  further  collisions  eventially 
spreading  their  energy  over  the  whole  nucleus.  For  a 
time,  the  nuclear  state  will  grow  increasingly  complex 
until  it  reaches  a  statistical  equilibrium.   This  state 
will  be  a  complicated  admixture  of  shell-model  and 
collective  configurations  involving  excitation  of  the 
many  degrees  of  freedom  of  the  nucleus.   Eventually 


enough  energy  will  be  concentrated  in  a  particle  or 
group  of  particles  to  allow  escape  from  the  nucleus 
and  upon  emission  the  nucleus  villi   decay  into  some 
new  excited  state.  This  process  is  referred  to  as 
scattering  via  the  formation  of  a  compound  nucleus 
and  is  usually  characterized  by  the  long  time  necessary 
for  its  completion,  which  is  about  a  thousand  times 
that  for  the  direct  interaction.   Since  these  two 
processes  exhibit  widely  different  cross  sections,  a 
method  for  separation  into  the  direct  interaction 
and  compound  nucleus  components  is  desirable.  No 
suitable  technique,  however,  has  been  found.   Instead 
reactions  are  chosen  which  are  predominently  either 
direct  interaction  or  compound  nucleus  types  and  where 
errors  made  in  neglecting  the  other  are  small. 

Direct  interaction  processes  involve  only  a  few 
of  the  many  degrees  of  freedom  of  the  nucleus.  The 
minimum  number,  in  fact,  would  be  just  those  necessary 
to  specify  the  ingoing  and  outgoing  channels  of  the 
reaction.  The  concept  of  the  direct  interaction  can  be 
extended  to  any  projectile  if  it  excites  only  one  degree 
of  freedom  in  the  nucleus   specifically  a  collective 
degree  of  freedom.  For  example,  it  would  certainly  be 
a  direct  interaction  if  the  projectile  bounced  off  a 
spheroidal  nucleus,  leaving  it  rotating.  These 
rotational  states  are  typically  low  lying  states. 


Compound  nucleus  reactions,  however,  involve 
many  degrees  of  freedom  due  to  the  complicated  nature 
of  the  nuclear  excited  state.   It  would  then  be 
fortunate  to  have  a  detector  capable  of  indicating 
the  number  of  degrees  of  freedom  in  the  wave  functions 
of  the  intercepted  particles,  but  alas,  none  has  ever 
been  found.   In  a  typical  compound  nucleus  reaction 
the  probability  of  a  particle  being  emitted  leaving 
the  nucleus  in  a  lovi   lying  state  is  greater  than  those 
leaving  the  nucleus  in  the  higher  excited  states. 
This  is  due  to  the  difficulty  of  barrier  penetration 
which  inhibits  the  emission  of  particles  with  energies 
just  above  the  separation  energy.   Figure  (1-1)  shows 
a  rough  plot  of  the  number  of  particles  emitted  vs. 
the  excitation  energy  of  the  nucleus,  where  €g   is  the 
ground  state  energy. 


Mill  iilil 


FIGURE   (1-1) 
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The  vertical  lines  show  the  discreet  spectrum.  The 
detector,  however,  does  not  possess  perfect  resolution 


and  the  energy  spectrum  it  sees  is  averaged  over  an 
interval  AS.   The  solid  curve  shows  such  an  averaging. 
Thus  it  is  seen  that  although  the  compound  nucleus 
would  rather  emit  particles  leaving  the  energy  of  the 
nucleus  near  its  ground  state,  the  majority  of  the  flux 
occurs  in  the  higher  lying  states. 

Hence  it  is  seen  that  inelastic  scattering  to 
low  lying  states  should  go  predominantly  by  the  direct 
interaction  process.   This  theory  will  be  developed 
in  the  following  chapters. 


CHAPTER  II 

THE  SCATTERING  CROSS  SECTION 

Introduction 

A  useful  concept  in  relating  theory  to  experiment 
is  the  differential  cross  section.   In  the  usual  picture 
of  scattering,  a  collimated  beam  of  particles  impinges 
upon  a  target.   The  number  of  particles  per  second,  dM, 
scattered  into  an  element  of  solid  angle  is  proportional 
to  the  incident  flux  N  and  the  element  of  solid  angle  cl-O.. 

dM  =  cr(e>4)  Wd-Sl-  (2"1) 

The  constant  of  proportionality  Cr(Q-,(j))  has  dimensions 
of  area  and  is  called  the  differential  cross  section. 
It  can  be  a  function  of  9  and  <j)  the  angular  coordinates 
of  d£L,   but  is  usually  only  a  function  of  9   taking  the 
polar  axis  as  the  beam  direction.  The  cross  section  is 
related  to  Schrodinger ■ s  equation  in  the  following  manner. 
At  distances  far  from  the  scattering  center,  the  wave 
function  describing  the  scattered  particles  can  be 
written  as  the  sum  of  a  plane  wave  corresponding  to  the 
incident  beam  and  an  out  going  spherical  wave  due  to  the 


scattering. 


-4-  =  k°®  ek'*'' 


(2-2) 


r 

The  function  i(&<&)   is  known  as  the  scattering  amplitude 
and  contains  all  of  the  angular  dependence  of  the 
scattered  beam. 

The  quantum  mechanical  expression  for  the  scattered 
flux  is  given  by 


*SCAr    S^ScAr       A.  scat  y  S^-SCAT 


(2-3) 


v/here  \mp  is  the  reduced  mass  of  the  scattered  system. 
Applying  this  to  (2-2),  an  expression  is  obtained 
giving  the  number  of  particles  being  scattered  /  unit 
area  /  sec. 

S^r    -  ^  p—  (2-4) 

v/here,  in  the  asymptotic  region  considered,  only  the 
radial  component  contributes  significantly  to  the  flux. 
The  number  of  scattered  particles  passing  through 
an  element  of  area  ciA  of  a  sphere  of  radius  r  located 
concentric  with  the  scattering  center  is  just  dU. 
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^   =(sScAr)dA  =  tkr   U(e$)l  dxi.   (2-5) 

mp  I      « 
where  d.CL   -   iA 

The  incident  flux  is  computed  using  an  expression 
similar  to  (2-3),  where  the  incident  wave  is  6. 

SlKit  =   ^JSt   39  Nl  (2-6) 

The  quantity  (In  is  then  divided  by  the  incident  flux 
and  the  resulting  expression  compared  with  (2-1)  to 
obtain  the  scattering  cross  section; 

a(G<t>)=  (-J(t)M2       (2-7) 

In  the  above  intuitive  argument,  it  should  be 
noted  that  any  interference  between  the  incident  and 
scattered  waves  has  been  ignored*   It  would  thus  appear 
that  (2-7)  is  inaccurate  and  of  limited  usefulness. 
The  incident  beam  however,  has  a  macroscopically  small 
cross  sectional  area  and  thus  would  only  interfere  with 
the  scattered  wave  at  both  very  small  and  very  large 
scattering  angles.  Since  the  density  of  particles  in 
these  beams  is  small,  they  will  not  greatly  influence 


one  another  even  in  these  cases.  For  a  more  rigorous 
treatment,  the  reader  is  directed  to  reference  £~^-J* 
where  it  will  be  noted  that  the  result  is  the  same* 
It  is  the  purpose  of  the  remaining  sections  of 
this  chapter  to  obtain  expressions  for  the  scattering 
amplitude,  using  the  Schrodinger  equation  with  an 
appropriate  potential.   As  this  is  a  study  of  the 
direct  interaction,  the  Hamiltonlan  will  exclude  all 
flux  that  would  go  into  populating  the  states  of  the 
compound  nucleus  and  would  itself  describe  only  those 
f ew  internal  degrees  of  freedom  pertinent  to  the  direct 
interaction.  Typically,  it  would  describe  the  motions 
of  the  one  or  two  nucleons  to  which  momentum  had  been 
transferred, 

i 
The  Schrodinger  Equation 

Consider  the  Shrodinger  equation  in  the  center 
of  mass  system.  The  subscript  f  means  "with  respect 
to  the  final  channel  f  •" 

—  — 

The  first  two  terms  on  the  left  are  the  kinetic  energy 
and  the  optical  interaction  potential  of  the  separating 
fragments,  and  ffC^,)   is  the  Hamilton! an  describing 
their  internal  motions.   VC&;.§r)  is  the  residual 
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interaction  energy,  with  respect  to  the  final  state 
channel,  which  gives  rise  to  transitions  to  excited 
states  of  the  residual  nucleus.  The  optical  potential 
is  of  the  form 


^~        (i  *  e(*-^) 


where  the  term  on  the  left  is  the  Coulomb  potential 
for  a  point  charge  2  and  an  incident  particle  of 
charge  ^   .   The  term  on  the  right  is  the  phenomenological 
complex  nuclear  potential  first  introduced  by  Saxon 
and  thus  bears  his  name.   Vi/hile  it  will  suffice  here 
to  give  a  plausibility  argument  for  the  complex  nature 
of  this  potential,  the  reader  is  referred  to  £~zj   and 
£~^J  for  a  theoretical  justification.  The  study  of 
dispersion  theory  in  physical  optics  has  shown  that 
absorption  of  light  in  a  cloudy  or  semi-transparent 
material  can  be  represented  mathematically  by  the 
imaginary  term  of  a  complex  index  of  refraction.   In 
the  study  of  nuclear  reactions,  the  phase  of  a  particle 
wave  impingent  on  a  potential  is  shifted  as  it  passes 
through  the  region  of  the  potential  and  it  is  known 
that  the  scattering  amplitude  is  completely  describable 
in  terms  of  these  phase  shifts.5  Completely  analogous 
to  the  optical  case,  if  these  phase  shifts  are  complex, 
the  imaginary  part  would  produce  absorption  of  part  of 
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the  Incident  beam.   It  is  then  the  imaginary  part 
of  the  complex  potential  which  gives  rise  to  absorption. 
This  term  literally  absorbs  out  all  flux  which  is  not 
described  by  the  potentialVG;?),  such  as  the  compound 
nucleus  flux,  thus  eliminating  it  from  the  scattered 
beam.   The  denominator  gives  the  radial  dependence  of 
the  potential  combining  its  short  range  nature  with  a 
diffuse  edge.   It  will  be  mentioned  in  passing  that 
many  variations  on  this  appear  in  the  literature. 
Hov/ever,  they  often  contain  so  many  parameters  that  the 
question  arises  as  to  whether  or  not  the  resultant  fits 
contain  any  physics  or  are  just  attempts  at  curve 
fitting.  The  variables  £  (or ^  in  the  incident  channel) 
are  those  few  internal  degrees  of  freedom  considered. 
The  Hamiltonian  7tf(  ^p)  is  assumed  to  have  eigen 
functions  V-(  £F)  which  are  solutions  of  the  equation 

To  solve  equation  (2-8),  It  is  first  written 
with  the  two  potential  terms  isolated. 

(-J£*+#&  "E)*  =  -Cm  wH        (2"10) 

To  remove  the  £_  dependence,  the  equation  is  multiplied 
by  the  complex  conjugate  of  V^-teVand  integrated 
over  the  internal  variables. 
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(2-11) 


Using  (2-9)  and  the  Hermitean  nature  of  /tvjfj* 
equation  (2-11)  can  be  written,  using  Dirac's 
"bra-ket"  notation  Tor  the  integrals  as 

fp  -   <V/p^  \  ^  -      (2-13) 

^-   ^(*-K>  (2-14) 

i  . 

Y^r  is  the  coefficient  in  an  expansion  of  *X  into  a 
set  of  eigen  functions  V^pO??)  an<i  can  be  thought  of 
as  a  projection  of  ^Ifon  the  basis  vector  for  the  final 

channel  "  f  ".   The  homogeneous  part  of  Equation  (2-12) 

Pi 
is  the  Helmholtz  wave  equation  with  solutions 

<p, -  e    — 
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General  Integral  Equation  Solution 

Equation  (2-12)  is  transformed  into  an  integral 
equation  through  the  use  of  a  Green's  function.   It 
can  be  shown  that 


%  =  -  2gr//f&^<ifr|U*v|4>^      ( 


2-15) 


is  a  solution  to  (2-12),  where  the  Green's  function 
KQCpjCp)  #  is  a  solution  to  the  equation 

(j7*+ **)  K(ev,t;)  =-£(rv-r=')     '  <««> 

as  long  as  r^^s  ry.   ,  (2-16)  is  just  the  wave  equation 
mentioned  earlier  with  outgoing  solution 


A  -  {J-  kpm  £> 


7 

This  is  expanded  in  a  series  of  Legendre  polynomials 


CO 


elk*  *Yl* fr+OftCk^fJs*^    (2"17> 


X*o 


7X 


Q 

Using'  the  Legendre  addition  theorem  ,  (2-17)  can  be 
expressed  in  the  space  fixed  coordinate  system. 
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&'**  4«T!JX 'u^X^X^   (2-i8) 


The  primed  and  unprimed  variables  are  the  angular 
coordinates  of  Jc^  and  j?^  respectively.  At  large 
distances  from  the  origin,  the  asymptotic  form 
of  j«(kr)  ma7  De  used. 


ft  (kpti  -*~  ^  Sl "  fe *  - &)     (2-i9) 


We  proceed  by  developing  a  series  solution  for  the 
Greeks  function.   All  "f "  subscripts  will  be  dropped 
until  they  again  become  necessary.  The  three 

dimensional  Dirao  delta  function  expressed  in  spherical 

9 
polar  coordinates  is 

The  spherical  harmonics  form  a  complete  set  over  the 
surface  of  a  sphere  and  a  mathematical  statement  of 
this  is  the  following  "closure  relation"' 
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The  total  delta  function  can  then  be  expressed  in1 terms 
of  this  series.  The  closure  relation  suggests  the 
following  expansion  for  the  Green's  function. 

This  is  now  inserted  into  the  defining  equation, 
(2-16),  where  the  Laplacian  in  spherical  coordinates  is 

Equation  (2-16)  becomes 


-a->  ^  .# 


The  spherical  harmonics  are  a  solution  to  the  f ollowing 
differential  equation. 
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This  is  solved  for  the  second  set  of  square  brackets 
and  the  two  series  are  equated  term  by  term. 


A'm^W)  (2-26) 


d 


The  radial  Green's  function  is  a  solution  to 

ififobjr  (g-MM)£(Cri=-±,Z(.Y--r')      (2-29) 

Consider  the  homogeneous  part  of  this  equation. 

if-P.(b-)  +  t£-4M)ttlkA  =0   .  (2.30) 


As  this  is  a  second  order  differential  equation, 
it  has  two  independent  solutions.   The  radial  Green's 
function,  Q?(^^)f   is  built  from  these  two  solutions 
by  the  following  prescription. 
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G(r,r') 


Region  I   •   Region  II 


FIGURE  (2-1) 


Figure  (2-1)  shows  a  plot  of  some  general  Gzifjt) 
as  a  function  of  r,   It  must  be  a  solution  to  (2-30) 
in  both  regions,  satisfying  the  boundary  condition 
at  the  origin  in  region  I  and  the  boundary  condition 
at  infinity  in  region  II,   While  it  is  required  that 
the  magnitude  of  Q»  C/;/-')  be  continuous  across  the 
boundary  at  rsr',  there  will  be  a  discontinuity  in 
the  slope.  The  magnitude  of  this  discontinuity  can  be 
obtained  by  integrating  equation  (2-29)  from  r-£  to  r'+e 
and  taking  the  limit  as  €-r0. 


iisA 


Limit 


<J(£$M) 


r-&. 


h'-t-G 


[If-  S/^qfcf'JJr 


4  £M 

or-*- 


* 


i&ii.  -  -# 


(2-31) 
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The  second  integral  contributes  nothing,  as  Qt(r,i'/) 
is  continuous  across  the  discontinuity  at  yl •  Two 
solutions  to  equation  (2-30)  will  be  obtained. 
Introduce  the  substitution 


fi(ki-)  -■f&Z(kl) 


This  is  inserted  into  the  homogeneous  equation  (2-50), 
yielding 


+  -|rS+O<-i*0L)?  =  O   Ca-32) 


This  is  Bessel's  equation  of  index  JL  +  k  and 


Z^t-)  =  TYlctO      . 

Multiply  through  by  the  constant  y/gr 

-C(b-)  m  kr  fip  JCIcO  =  Icr-t  Cfehi    (2-53) 

where  i#(*^*")  is  the  spherical  Bessel  function. 
Consideration  of  equation  (2-17)  shows  this  to  be 
kr  times  the  radial  wave  function.  Here,  ^(kr)  is 
regular  at  the  origin  and  thus  satisfies  the  boundary 
condition  in  region  I.  The  following  relation  exists 
between  the  spherical  Bessel  and  Hankel  functions.1 

(A.     A 

(2-54) 


pU^-.±-($M+\Uty 
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The  functions  ft*    and  ft*    are,  respectively,  tne 
spherical  Hankel  functions  of  type  one  and  type  two. 
They  are  both  solutions  to  equation  (2-50). 

id 

Asymptotically  ffi    goes  as  an  outgoing  spherical  wave 

j(Z) 

and  fig    as  an  incoming  one.   As  r  approaches  infinity, 

there  can  be  only  outgoing  scattered  waves.   Hence 
a  solution  satisfying  the  boundary  condition  in 
region  II  is  s?Akr)  •  The  Green's  function  is  written 

C,  trailer)        r<i~l    • 

sfc^H      >  C2-35) 


% 


i[kr^V)]-Qi^,^) 


(2- 31a) 


subject  to  the  condition  imposed  by  (2-31) . 

0) 

The  Wronskian  relation  for  these  two  functions  can  be 
used  to  evaluate  the  constants  C]_  and  Cg. 
Sturm- Li ouvi lie  theory12  gives  the  Wronskian  of  two 
independent  solutions  to  equation  (2-30)  as 

>0), 


Wfc 


>r),g^))  =  a  <2-36) 


where  A  is  a  constant . 
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As  this  is  true  for  all  r,  the  asymptotic  forms 
of  these  expressions  may  be  used  in  computation. 
The  asymptotic  form  of  is/S^^   was  already  given 
in  (2-19)  and  for  large  r 

Calculation  of'  the  V/ronskian  yields 

A  =  k 


A  comparison  of  (2-31a)  with  (2-37)  shows 

0). 


C,   =  k>r") 


Cz  = 


=  fa^ 


The  Green's  function  becomes 

Y^    is  the  smaller  of  r  and  r' 
f^  is  the  greater  of  r  and  r' 


(2-39) 


21 


Combining  equations  (2-23),  (2-26),  (2-28),  and 
(2-39),  the  total  Green's  function  is  obtained. 

^„  m-4  •  (2- 23a) 

Equation  (2-23a)  is  then  used  in  (2-15)  to  give 
the  solution  for  ip. 

The  General  Scattering  Amplitude 

As  was  discussed  in  the  introduction  of  this 
chapter,  the  scattering  amplitude  is  the  coefficient 
of  the  outgoing  scattered  spherical  wave.   Equation 
(2-15)  will  be  used  to  find  the  asymptotic  form 
of  Tp-.   Equation  (2-23a)  can  be  simplified  for  large  r, 
r  being  )y  r,  where  K  can  be  replaced  by  its 
asymptotic  form  given  in  (2-39).  The  Green's  function 
becomes 

where  the  primed  angular  coordinates  refer  to  K.   and 
the  unprimed  ones  to  j^-v/hich  is  taken  to  be  along 
the  vector  kp.   Consider  the  expression  in  (2-18). 
Taking  the  complex  conjugate  yields 


22 


s> 


where  the  primed  coordinates  refer  to  J^.' and  the 
unprimed  ones  to  JTp. 

r 
Using  equation  (2-15) 

The  plane  wave  <Q  ~-~-±3   a  solution  to  the  homogeneous 
part  of  (2-12)  and  can  always  be  added  to  the  scattered 
solution.  The  scattering  amplitude  is  then  equated 
with  the  coefficient  of  -=■  G 

ft^)  =  -^XvF^^lu+vl±^  (2"40) 

where  the  Dirac  notation  is  again  used  to  indicate 
integration  over  all  internal  and  external  variables. 


The  Elastic  Scattering  Amplitude 

In  a  description  of  only  the  elastic  component 
of  the  scattering,  the  potential  T/C£  %.)   which 
describes  transitions  to  excited  states  is  set  equal 
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to  zero  and  the  interpretation  or  the  imaginary  part 
of  the  complex  potential  is  broadened  to  include 
absorption  of  the  entire  inelastic  flux.  The  initial 
and  final  channels  are  then  the  same  and  the  outgoing 
wave  function  ^  is 

T=  *<&t%£  (2-41) 

where  Vr(^r)  and  JM  ^ryJ3r)  a^e*  respectively,  the 
internal  and  elastic  scattering  wave  functions  for 
the  initial  channel*  £,  satisfies  the  equation 

The  elastic  scattering  amplitude  is  obtained  from 
(2-40) 


Inelastic  Scattering  -  Distorted  V/aves 

The  direct  interaction  inelastic  scattering 
component  is  to  be  isolated  from  the  scattered  wave 
and.  using  techniques  developed  in  the  first  part  of 
this  chapter,  the  inelastic  scattering  amplitude  v/ill 
be  found.   Equation  (2-8)  is  re- written  with  the 
potential  which  describes  the  inelastic  transition 
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isolated. 


Using  (2-9)  and  the  Hermitean  nature  or  ^  (  £  ), 
the  ^.dependence  is  removed  in  the  same  manner  as 
for  (2-12) 

(y2-2^£H  =  ^«t|Vft^)|t>   (3-43) 

2^  =  ^UCrP)  (2-44) 


Kp  and  Tp  are  discussed  after  equation  (2-12)  and  v/ill 
be  reiterated  here  for  the  readers'  convenience. 

The  homogeneous  part  of  (2-43)  has  as  one  of  its 
eigen  solutions,  the  elastic  scattering  wave  function 
discussed  in  relation  to  equation  (2-41).  A  series 
solution  corresponding  to  an  outgoing  wave  will  be 
obtained.   Consider  the  homogeneous  equation  with 
the  f  subscript  dropped  as  before, 

'    (372-2^+  £)X%  O         '    (2-45) 
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The  form  of  the  desired  solution  is  such  as  to  reduce 
to  expression  (2-18)  when  the  optical  potential  is 
removed,  as  certainly  (2-45)  reduces  to  the  Helmholtz 
wave  equation.  As  the  optical  potential  is  a  function 
of  only  the  radial  variable,  JC  (J£,JT)  will  have 
azimuthal  symmetry  if  Is   is  assumed  to  be  along  the 
z-axis.   (Note  expression  (2-17)).  The  discussion 
following  Equation  (2-8)  gave  the  optical  potential 
to  be  of  the  form 

A  plausibility  argument  for  the  complex  form  of  the 
nuclear  potential  was  given  as  well  as  an  expression, 
typical  of  those  appearing  in  the  literature.  No 
explicit  form  of  Vjj  is  necessary  in  the  following 
calculations.  However,  it  will  be  assumed  that  the 
nuclear  forces  are  short  ranged  and  thus  ,Vjj  will 
rapidly  approach  zero  for  distances  greater  than  the 
nuclear  radius.  Equation  (2-45)  becomes 

[y*  +  |*  -4*4  _  ^  Usl]^  O   (2-45.) 
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Coulomb  Wave  Functions 

Consider  the  special  case  of  zero  nuclear 
potential.   This  situation  would  be  encountered  if 
the  impingent  beam  lacked  sufficient  energy  to 
overcome  the  coulomb  barrier  and  penetrate  into  the 
nuclear  material.   The  eigen  solutions  become  coulomb 
wave  functions,  satisfying  the  equation 

[V 2  +  £  -  ^J?]  §££)   =  O  (2-47) 

The  second  chapter  of  reference  £"?>J   is  an  excellent 
source  for  information  regarding  solutions  of  this 
equation.   It  turns  out  that  (2-47)  has  a  regular 
solution  of  the  form 

•   (bc  =-^"  Z-fO-2)  (2-48) 

If  z  is  the  polar  axis,  then  z  =  r  Cos9  and  f  (r-z)  has 
azimuthal  symmetry.  Using  (2-24),  the  Laplacian  is 
expanded  in  spherical  polar  coordinates  and  the 
following  substitution  is  made. 

U=  V-O-Cos©) 
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The  substitution  is  carried  out  using  the  chain  rule. 

ox  du  3x 

The  results  in  the  following  differential  equation 
for  fCU) 

U  ^.     +(l-UcU.j jL- Dclcjftu)- O   C8.50) 

The  factor  k  is  put  completely  within  the  variable 
with  the  transformation 


-j*  Cku. 


(2-51) 


Equation   (2-50)   then  becomes 

V4L+   0-V)ft,    +  <*"]&)-  G  (2-52) 


[' 


Equation  (2-52)  is  in  the  form  of  the  confluent 
hypergeometric  equation  which  is  solved  in  appendix  A, 
A  comparison  of  equation  (A-l)  and  its  solution 
(A-4)  with  (2-52)  gives 


-fa;  =  x  (-kxiiiv; 


and 
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(J>   =  Ce  ,  Fi  ("Cct  1 1 1  I  lcO-2))        (2-53) 


A  consideration  of  equation  (A-4)  shows  it  to  be 
non- singular  at  the  origin  and  therefore  the  solution 
6C  is  the  regular  coulomb  wave  function.  The  constant 
C  is  determined  so  that  asymptotically 


"i  -  e     +  * 


SCAT- 


where  as  usual  6  scat  is  the  scattered  outgoing 
spherical  wave.   In  appendix  A,  the  following  asymptotic 
solution  (A- 19)  is  derived. 


*uw~  &**%&& 


.a.-L 


making  the  substitutions,  equation  (2-53)  becomes 

to  Cm)  rO-tcO 

the  logarithm  is  a  slowly  varying  function  with  respect 
to  its  argument.  Hence  at  large  distances  from  the 
scatterer,  the  first  of  these  two  terms  is  essentially 


29 


an  outgoing  plan©  wave  and  the  second  an  outgoing 
spherical  wave.  The  coefficient  of  the  plane  wave  term 
was  defined  as  being  unity  and  thus  the  value  of  C 
is  determined. 

Equation  (2-47)  is  also  separable  in  spherical 
coordinates,  thus  leading  to  a  decomposition  of  A) 
into  partial  waves, 

Let  rf>   =.   FCP)  (&&) 

When  this  is  substituted  in  (2-53)  two  equations  for 
P  and  ©  are  obtained. 

S»3  i<^  $D  +Jia+>)  ®  -  °       (2-5S) 

where  ^(^C-hl)  is  the  separation  constant.  Equation 
(2-56)  is  Legendre!s  equation  which  has  Legendre 
polynomials  as  its  solutions. 

<B>  =  ^  (Cos©) 

In  equation  (2-51)  let 

p     =     -  ^c  fcr        rf         (2-58) 
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V/ith  this  transformation,  (2-57)  becomes 

Ut     §■  =*(%)***  Vrf)  (2-60) 

After  some  manipulation,  the  equation  for  Vo.  becomes 

f> &  +  (ai+a^M  _  tf+  /  +u)  l£  =  o      (2-61) 

A  comparison  of  (2-61)  to  (A-l)  shows  Vn    to  be  a 
solution  of  the  confluent  hypergeometric  equation 

where , 

\jl  a, If  (jUl+u  1  aX+al^)  (2.62) 


The  radial  coulomb  wave  function  becomes 

(■[&*■)  =  /^e  ?fcf)£+', /■"&■/  +  c-<x|22+2|-56i:f-)   (2-63) 

where  the  Ao  are  determined  by  the  same  asymptotic 
considerations  as  before,  namely  that  d>c  behaves  as  an 
outgoing  plane  wave  at  large  r.  The  asymptotic 
expression  for  F?  is  calculated  in  the  same  manner 
as  for  equation  (2-53).  By  substitution  of  the 
proper  parameters  into  (A- 19),  the  following  result 
is  obtained. 
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(2-64) 


P(£-H-tCc^) 

The  general  theory  of  complex  variables  gives 

where  .Cjf  a.  arg(r  (Jj+l -K.0O ")  (2-65) 


with  this  substitution,  expression  (2-64)  simplifies  to 

A  consideration  of  the  plane  wave  expansion 
in  (2-17)  shows  the  corresponding  radial  wave  function 
to  be  jfl(kr),  with  its  asymptotic  form  given 
in  (2-19).   The  constants  kn   are  chosen  so  that  these 
two  functions  approach  each  other  as  \r  — *■ 


The  total  solution  for  <j)c  is  then  the  sum  of  all  the 
partial  waves.  This  is  shown  by  using  expressions 
(2-53)  and  (2-55)  for  <J)C  and  the  contour  integral 
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representation  of  the  confluent  hyper geometric 
function  given  in  (A-7) 

-  o-e^)  p(-u) ' '  *  e      *    **  **  C2'68) 

using  equation  (2-17),  the  exponential  e       'is 
expanded  into  a  series  of  Legendre  polynomials 
where  z  =  r  Cos© 


Cx> 


<£-<-  -  2_l*fa+t)  \(kh)  %  (Coy^j  (2-69) 

Where 

Kw*(ijr)mp  $***¥**?*     {2-70) 

a  consideration  of  equations  (2-30),  (2-33),  (2-57), 
(2-63),  and  (2-67)  with   equal  to  zero  in  the  last 
three,  gives 

^^J»dMfVlF(wNh2u)   (2-71) 

using  the  series  form  of  the  confluent  hypergeometric 
function  (A-4),  equation  (2—71)  becomes 


* 


(2-72) 


Letting  x  =kr(l-t),  (2-72)  is  substituted  into 
(2-70)  and  the  order  of  summation  and  integration 
is  reversed.  This  gives 

CJLi— 


V*   (2-73) 


*£«^^wr 


The  integral  given  in  the  parenthesis  is  a  contour 
integral  representation  of  the  beta  function  as  given 
in  expression  (A-6) 

with  this  substitution,  (2-73)  becomes 

0  ^ 

By  equation  (A-4),  this  is  written 
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(2-74)  is  compared  with  (2-63)  and  (2-67),  giving 

£(fcr)  m,  %(&)'■£  (2-75) 

where  C£  is  defined  in  expression  (2-65).   Expressions 
(2-69)  and  (2-75)  are  combined  and  the  partial  wave 
expansion  of  (j)  is  obtained 


Optical  Wave  Functions 

The  eigen  solution  to  equation  (2-45a)  corresponding 
to  an  outgoing  wave  will  now  be  obtained.   It  was  shown 
in  the  previous  section  on  coulomb  waves,  that  if  the 
nuclear  potential  is  neglected,  the  equation  is 
separable  in  spherical  polar  coordinates,  and  since  V^ 
is  a  function  of  only  the  r  variable,  the  separability 
is  unaffected.   If  the  Laplacian  is  expanded  and  the 
variables  separated,  the  angular  wave  function  is 
unchanged.  The  equation  for  the  radial  coordinate  is 


35 


(£x+   £-^-^-M0)^).O    (2-77) 


The  total  solution  is 


%*(kr)-m  -l^.  (*$(**) %&**)         (2-78) 


where,  as  usual,  the  constants  C^  are  determined 
from  the  asymptotic  considerations. 

It  was  shown  in  equation  (2-64)  that  the  regular 
asymptotic  form  of  the  radial  coulomb  wave  function 
could  be  built  from  a  linear  combination  of  the  two 
asymptotic  exponential  solutions 

e  x 

where  the  plus  and  minus  signs  correspond  to  outgoing 
and  ingoing  waves,  respectively •  As  the  nuclear 
potential  was  postulated  as  being  short  ranged,  these 
may  also  be  used  in  determining  the  asymptotic  form 
of  ^(kr) .  This  is  seen  by  substituting  the 
exponentials  into  equation  (2-77)  with  the  following 
result • 
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It  is  seen  that  as  long  as  V-^   approaches  zero  at  least 
as  fast  as  v     ,   an  asymptotic- form  of  j?^(kr)  can  also 
be  formed  from  these  exponentials  with  error  certainly 
no  greater  than  that  inherent  in  the  coulomb  solutions. 
The  phenomeno logical  nuclear  potential  postulated  in  the 
earlier  part  of  this  chapter  meets  this  requirement* 
Prom  consideration  of  the  corresponding  coulomb  solution, 
the  regular  asymptotic  solution  to  the  radial  optical 
wave  function  is 

jfa&ir)-**  Sin0tta-ef4h3lcir-4flr -*$r-* 4.)  (2-79) 

where  <?£.is  an  additional  phase  shift  which  the  wave 
would  experience  in  passing  through  the  region  of  the 
nuclear  potential.  The  constants  Cji  are  now  determined 
from  the  following  two  considerations:   (1)  The  nuclear 
potential,  being  short  ranged  cannot  affect  the  incoming 
asymptotic  solution  and,  (2)  When  %  is  set  equal  to 
zero  *)£?*  must  reduce  to  the  coulomb  wave.  As  it  was 
shown  that  the  solution  could  be  written  as  a  linear 
combination  of  the  ingoing  and  outgoing  coulomb  waves, 
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the  first  requirement  becomes 

^  ^  (2-80) 

The  solution  of  this  equation  yields 


Ae.  --  ai  4 


L<U 


C*  =  C&  Q 


I  Jo> 


(2-82) 


In  order  for  (2-78)  to  reduce  to  equation  (2-76) 
when  Vjj  and  thus  oN  are  zero, 

e>    £*fci*rj«4*  (2.83) 


and  the  total  wave  function  becomes 

if^&T^^iwfa1^      (2-84> 

where         CP-nj.  =  d^-K^fi.  (2-85) 
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Using  the  Legendre  addition  theorem  to  write  this 
in  terms  or  the  angular  coordinates  of  the  vectors 
j^  and  r,  equation  (2-84)  becomes 


8.**     K>»—JL 

<£-  (2-S6) 


where  the  primed  variables  will  be  taken  as  the 
coordinate  corresponding  to  k. 

Inelastic  Integral  Equation 

The  following  development  parallels  that  used 
in  the  study  of  the  general  integral  equation,  the  main 
difference  being  that  distorted  waves  are  used  instead 
of  plane  waves.  The  inelastic  component  of  the 
scattering  has  already  been  isolated  in  equation  (2-43) 
and  the  corresponding  integral  equation  for   is 

jyY(r,r!)»  the  Green!s  function  is  a  solution  to.  the 
following  equation.  Note  that  the  f  subscripts  have 
been  dropped  as  before. 
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It  was  seen  in  equations  (2-21)  and  (2-22) 
that  the  Dirac  delta  function  could  be  expressed 
as  a  series  of  spherical  harmonics  as  follows 

Equation  (2-88)  differs  from  (2-16)  in  the  previous 
treatment  only  in  the  addition  of  a  potential  term 
which  contains  no  angular  dependence.  The  series 
expansion  of  the  Green fs  function  defined  by  (2-88) 
will  then  have  the  same  angular  dependence  as  (2-23). 
and  can  be  written 

'(A/hen  this  is  inserted  into  (2-88)  using  the  Laplacian 
given  in  (2-24),  the  radial  Green* s  function  is  seen 
to  be  a  solution  of  the  equation 

, .^hfc)^-«rt-^)Mt-  -£ Jow4  (2.90) 

The  reader  is  referred  to  equation  (2-25)  and  its 
subsequent  discussion,  where  more  detail  is  given 
to  a  similar  calculation. 

Letting    gCt^r)     -     h-^CcH)  (2-9!) 
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equation   (2-90)   becomes 

£c£  +  £  _•*  Lr)  -  MM  \%(r>r>)  =  -XfcrwO    (2-92) 

The  radial  Green's  function  is  obtained  from  two 
independent  solutions  to  the  homogeneous  part  of  this 
equation,  one  regular  at  the  origin,  and  the  other 
corresponding  to  a  purely  outgoing  scattered  wave 
satisfying  the  boundary  condition  at  infinity.   The 
regular  solution J?fct) is  obtained  from  equation  (2-77) 
v/ith  its  asymptotic  form  given  in  (2-79). 

Let  j^^be  the  solution  to  (2-77)  and  thus  to  the 
homogeneous  part  of  (2-92)  with  outgoing  boundary 
conditions.   It  was  shown  in  equations  (2-80)  and  (2-81) 
that  the  asymptotic  regular  solutions  could  be  written 
in  terms  of  the  two  irregular  exponential  solutions 

e 

where  the  plus  and  minus  signs  correspond,  respectively, 
to  outgoing  and  incoming  waves. 
Thus,  asymptotically, J^(k,r)  becomes 

titMr-^^-^+oh]         (2_93) 

K(kf)  -r  -a  (2  93) 
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Tiie   Green's  function  is  written 


%W  =  i 


6,  3i(k}i~)        r<<-' 
€zHji(kh)        r?h' 


(2-04) 


where,  analogous  to  equation  (2-31a),  the  discontinuity 
in  the  derivative  of ^{r, r')  evaluated  at  r'  is  written 

^L^Ckh}  -  *.$&*»  =  ~k     (2-95) 

The  constants  G--^   and  £g  are  determined  from  the 
Wronskian  relation 


~W(.&X)*-A' 


(2-96) 


where  A1  is  a  constant.  The  Wronskian  is  evaluated 
using  the  asymptotic  forms  of  J^  and  yj£  and  found 
to  be 


%4zJ&  -  mAjz  =  *t 


^w^ 


%   drv* 


(2-97) 


A  comparison  of  equation  (2-97)  with  (2-95)  yields 
the  values  of  the  constants. 

I 


61-  ir,-^^ 


&2    kl~  *■ 


(2-98) 
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Using  the  notation  of  (2-39),  the  radial  Green's 
function  becomes  ^ 


fQrt  *   i  CjfCkGl%&$)  (2-94a) 


and  the  total  Green's  function  is  obtained  from 
(2-89),  (2-91)  and  (2-94a). 

X  YU'4)  Y(?4)' 

This  is  then  used  in  expression  (2-87)  to  obtain  the 
outgoing  scattered  wave. 

Inelastic  Scattering  Amplitude 

Using  (2-87),  an  expression  f  or  **£ in  the  asymptotic 
region  is  sought.  For  large  r,  r>>r» ,  JfyfcfJ can  be 
replaced  by  its  asymptotic  form  given  in  (2-93)  and 
the  Green's  function  written. 

c  *f<ttr  (2-99) 

where  OV-j^is  defined  in  expression  (2-85)  and  r^  is 
taken  along  the  direction  of  kf  •> . 

Consider  the  following  part  of  (2-99) 
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Upon  comparison  with  (2-8S),  these  are  recognized 
as  being  similar. 

Expression  (2-100)  is  identified  with  the  complex 
conjugate  of  the  ;"time  reversed"  optical  wave  function, 
and  is  denoted  by^L.(kr).   In  the  "time  reversed" 
scattering  problem,  the  initial  and  final  channels 
have  exchanged  roles  and)£6<r)is  a  solution  to 
equation  (2-45a)  with  incoming  boundary  conditions, 
where  V*j  has  been  replaced  by  its  complex  conjugate 
to  simulate  flux  re- inserted  into  the  beam.   The 
function  is  perhaps  best  generated  using  the  previous 
calculations  for^caf:'(k,r)  and  the  Wigner  Relation 

f%ti^-t\-±£)  (2-101) 


Using  equation  (2-86), 

.-if-  ■ 


14 
The  parity  of  the  spherical  harmonic  is 


(2-102) 
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Hence,  equation  (2-102)  becomes 

which  agrees  with  (2-100) 

The  Green's  function,  (2-99),  is  then  re-written 
incorporating. the  "time  reversed"  wave  function. 

•   j&>*  -44e        ^  ^  (2"99a) 

Using  equations  (2-87)  and  (2-99a),  the  asymptotic 
expression  for  N^  is  obtained. 


YF  -  }J;^  -  e         ^^jm>£  (2-104) 

The  first  term  is  the  elastic  scattering  or  optical 
wave  function  which  is  a  solution  to  the  homogeneous 
equation  (2-45)  and  is  necessary  should  the  initial 
and  final  channels  coincide.  This  would  be  the  case 
of  elastic  scattering.  The  inelastic  scattering 
amplitude  is  the  coefficient  of  the  outgoing  spherical 
wave  • 

ty* = -»&<^iyi?b>    (2-ios) 

With  the  inelastic  component  of  the  scattering 
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amplitude  separated,  the  total  scattering  amplitude 
can  be  written,  not  only  as  in  equation  (2-40), 
but  also  as  a  sum  of  its  elastic  and  inelastic 
components.  Using  equations  (2-42)  and  (2-105), 
the  total  scattering  amplitude  is  written, 


t-kf-iri , ,  i . ,  vw\ 


(2-106) 


Matrix  Element  Formalism 

Customarily,  the  scattering  amplitude  Is  not 
used  In  the  literature,  but  rather  the  matrix  element 
giving  the  probability  of  transition  from  an  initial 
to  some  final  scattering  state,   A  consideration  of  the 
expressions  for  the  scattering  amplitude  shows  the 
matrix  element  to  be  given  by 


where  the  subscripts  i  and  f  refer  to  the  initial 
and  final  scattering  states,  respectively. 

In  this  formalism,  equation  (2-106)  becomes 

15 
the  Gell-Mann-Goldberger  relation  . 

which,  while  derived  here  for  the  specific  case 
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of  the  optical  wave  functions,  is  quite  general  and  is 
usually  obtained  from  a  unitary  transformation  of  the 
general  scattering  amplitude  matrix  element, 

TFS=  <vFe     ^lu+vi^)  (2-io9) 

v/here  IX  and  V  are  any  two  potentials  that  have  been 
isolated. 

The  angular  distribution  or  cross  section  is 
written  in  this  formalism  using  equation  (2-7) 

AV 

The  summation  sign  indicates  a  twofold  operation; 
(1)  due  to  the  inability  of  the  detector  to  distinguish 
betv/een  spin  orientations,  the  projections  of  the 
final  state  spins  are  summed,  and  (2)  the  projections 
of  the  initial  state  spins  are  averaged  as  when  time 
and  thus  the  direction  of  particle  motion  is  reversed, 
the  direction  of  any  spin  must  also  be  reversed  to 
preserve  time  reversal  invariance.    It  is  seen  that 
if  these  projections  are  averaged,  producing' an  "average 
spin"  initial  state  wave  function,  the  problem  is 
averted.   A  suitable  modification  would  have  to  be 
made,  however,  if  polarization  experiments  are  to  be 
considered. 
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It  now  only  remains  to  specify  a  suitable  potential 
and  scattering  wave  function  and,  using  (2-110),  the 
cross  section  can  be  obtained.  This  problem  is 
considered  in  the  following  chapter. 
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CHAPTER  III 
EVALUATION  OP  THE  DIFFERENTIAL  CROSS- SECTION 
FOR  COLLECTIVE  EXCITATIONS 
Introduction 

An  exact  calculation  of  the  matrix  element  T^.. 
is  at  best,  a  difficult  problem.   Transforming 
Schrodinger » s  equation  to  its  integral  equation 
counterpart  has  not  given  its  solution,  but  has  only 
provided  a  convenient  way  to  obtain  the  scattering 
cross  section  once  the  wave  function  has  been  found. 
The  specific  interaction  potential  is,  in  general, 
not  known  and  thus  an  exact  solution  for  x  is  not 
easily  found.  For  example,  if  the  angular  momentum 
transferred  in  the  reaction  is  non-zero,  the  interaction 
potential  will  not,  in  general,  be  a  scalar  function 
of  position  and  normal  methods  for  solution  of  the 
Schrodinger  equation  fail.  Physically  meaningful 
results,  however,  may  be  obtained  v/ith  appropriate 
approximations  for  the  wave  functions  and  the 
interaction  potential.   It  is  the  purpose  of  this 
chapter  to  discuss  the  techniques  involved,  and  to 
show  their  applicability  by  fitting  experimental 
angular  distributions* 

In  each  of  the  expressions  for  Tfi,  four 
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' integration  variables  are  employed.   They  are 
£.,  £j.f     ^*,    and  ^  where  each  is  a  vector  quantity. 
Undoubtedly,  both  the  initial  and  final  position  ■ 
variables,  and  the  initial  and  final  internal  variables 
are  related,  but  the  exact  nature  of  this  relation 
is  unknown.   It  would  be  convenient  to  have  but  one 
position  variable  and  one  internal  variable  when 
performing  these  integrations,  hence  for  want  of  this 
exact  relation  let 

*  (3-D 

*?x  -  £P  *   ? 

This  is  called  the  zero-range  approximation. 

The  angular  distribution  is  now  evaluated  using 
the  general  scattering  matrix  element. 

TFx=  <vPe~p"r|U+vi^._>  (s-2) 

It  will  be  assumed  that  the  wave  function  can  be 
at  least  partially  factored  and  takes  on  the  form 

X+) 


where  the  wave  function  v- ( £ )  is  the  eigen  function 

ft) 


for  the  initial  nuclear  state,  and  /^U^S.*?)  is  a 
solution  to  the  equation 
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It  should  be  noted  that  the  Hamilton! an  describing 
the  internal  motions  is  not  present.   Hence  the  & 
variables  can  be  treated  as  constants  and  the  function 

fz      depends  on  them  only  in  a  parametric  sense. 
T~.  is  then  obtained  by  a  two  step  calculation; 
(1)  equation  (3-3),  a  function  of  only  the  r  variable, 
is  solved  and  the  resulting  solution  used  to  compute 
the  transition  amplitude  for  a  fixed  value  of  *J , 

where   Tpx<$  =r  C£^\  tUV|  J^}  <3-5) 

and  (2)  this  result  is  multiplied  by  v1(^)  and  the 
projection  of  their  product  on  the  final  state  f  is 
found. 

Physically,  the  above  corresponds  to  the  assumption 
that  the  initial  nuclear  state  described  by  the  function 
v«(2s)  changes  so  slowly  with  respect  to  the  transit 
time  of  the  projectile,  that  the  scattering  wave  function 
is  completely  determined  by  the  value  of  v^(^)    at  the 
instant  of  the  reaction.  This  is  analogous  to  "stop 
motion"  photography,  where  a  moving  object  has  its 
motion  frozen  on  film  through  the  use  of  a  short  duration 
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flash  of  light.   In  the  time  considered,  the  fiT 
variables  are  almost  stationary  and  thus  it  is  seen 
that  tho  Hamiltohian,  ^(fe7),  can  be  neglected  in.  the 
integral  equation  leading  to  (3-2)  and  in  (3-3), 
These  assumptions  are  most  valid  when  high  energy 
projectiles  are  used  and  when  the  nuclear  excited 
states  are  collective  states.   Since  the  excitation 
energy  must  be  shared  with  the  fairly  large  mass 
of  the  nucleus,  the  collective  variables  exhibit  a 
slow  change. 

The  adiabatic  theory  (Diffraction  Scattering) 
will  now  be  developed  and,  using  the  computed  angular 
distributions,  will  be  compared  with  experimental 
alpha  particle  scattering  data. 

Diffraction  Scattering 
The  Praunhofer  Approximation 

The  calculation  is  most  conveniently  made  in  the 
scattering  amplitude  formalism;  hence  (3-5)  becomes 


J}fL 


3  -CJSr-t/ 


Equation  (3-3)  is  multiplied  through  by  2m/h  and  the 
resulting  expression  is  substituted  into  (3-6)  thus 
transforming  away  the  dependence  on  the  interaction 
potential. 
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f(W>*^ 


(*) 


(3-7) 


wnere 


This  is  put  into  a  more  convenient  form  by  recognizing 

— Ak  •  r 
that  the  exponential  e  -f'-is  a  solution  to  the 

wave  equation. 

V  e        =wp-e  - 


If  each  side  is  multiplied  by  ^£^and  the  right  and 
left  sides  are,  respectively,  added  to  and  subtracted 
from  the  integrand  of  (3-7),  the  result  will  leave  the 
integral  unchanged,  but  allow  it  to  be  written  in 
the  form 

Making  the  adiabatic  approximation,  k^kf,  the  first 
of  these  integrals  becomes  vanishingly  small  and  will 
be  neglected  in  the  remainder  of  this  development. 
Asymptotically  7%  can  be  written 
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where   'y  scat  *s  tilat  component  of  the  wave  due  to 
the  scattering.   This  is  inserted  into  the  second 
integral  or  (3-8)  v/ith  the  result 


(3-10) 


The  first  of  these  integrals  vanishes  as  in  (3-8)  when 
the  adiabatic  approximation  is  applied.   Green's 
theorem  relating  the  volume  and  surface  integrals 
of  potential  functions  is 

If  it  is  assumed  that  the  potential  TJ  -{"  V  is  non- 
vanishing,  only  within  some  finite  volume,  then  the 
second  integral  in  (3-10)  is  transformed  into  the 
following  surface  integral 

S" 

where  the  surface  S  is  chosen  so  that  the  interaction 
potential  vanishes  outside. 

The  impingent  particles  are  assumed  to  be  fairly 
energetic  so  that  the  high  energy  approximation  can 
be  used.  In  this  approximation,  the  paths  of  the 
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incoming  particles  are  assumed  to  follow  straight 
lines  parallel  to  the  vector  k. .   The  surface  S 
is  then  conveniently  chosen  to  be  a  cylinder  whose 
axis  lies  along  the  vector  k*  as  shown  in  figure  (3-1). 


O 
5 


■*■* 


FIGURE  (3-1) 

Y/ith  this  choice  of  surface,   "^scat  vanishes 
everywhere  except  in  the  shadow  plane  in  the  region 
of  the  nuclear  shadow.  The  entire  contribution  to 
the  integral  in  (3-12)  comes  from  this  same  region 
and  the  element  of  area  d©"  becomes 

d<r*  kT  dA 

mm 

\ 

where  k^  is  a  unit  vector  along  k^  and  dA  an  element 

of  area  in  the  shadow  plane.  Equation  (3-12)  is 
then  written 


Shadow 
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To  compute  the  Fraunhofer  amplitude,  the  nucleus 

is  assumed  perfectly  absorbing  or  "black"  to  the 

incident  beam  and  thus  "^2:  ^  0  since  there  can  be 

no  particles  in  the  shadow  region.  Oh  Is 

/  s  c  a  x>      , 

obtained  from  (3-9)  by  setting  Ife-  0. 

jscat   -  -  cr 

This  is  incorporated  in  equation  (3-11)  which  becomes, 


f^~  kc^p-^ 


(3-14) 


where  ©,  the  scattering  angle,  is  the  angle  between 
the  vectors  k*  and  Jo. .   The  vector  difference  kj_~kf 
is  the  momentum  transferred  to  the  nucleus.   It  has 
magnitude  2k  Sin  (  &/2)    if  the  adiabatic  approximation 
kr55  kf  *  ^  is  ma<ie  anci  i'or  small  scattering  angles 
has  a  direction  very  nearly  perpendicular  to  k^. 
Thus  for  small  ©,  the  momentum  transfer  vector  lies  »: 
in  the  shadow  plane  and  equation  (3-14)  can  be 
written 

r         .  /     v  ff  *fSnfiG»<(> 

TP^=^Cl+<is^)//e     *  hcWf>    (3-15) 

v/here  the  elemental  area  dA  is  written  in  plane 
polar  coordinates,  and  q>  is  the  angle  between  the 
momentum  transfer  vector  and  r# 
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The  Theory  of  Blair  and  Drozdov2'3 

The  collective  variables  ^r  can  only  affect  the 
scattering  amplitude  by  changing  the  size  and  shape 
of  the  nuclear  surface  and  thus  the  nuclear  shadow. 

Hence  the  collective  model  proposed  by  Bohr  and 

4 
Mottelson  is  appropriate. 

It  is  assumed  that  the  nuclear  surface  can  be 

described  by  the  following  multipole  expansion 

r=  Po+rw^  (3-is) 

V-H 

» 

where  the   *?tm  fl-re  the  deformation  distances  of  the 
LM  multipoles.   If  the  incident  beam  direction  lies- 
along  the  z-axis,  the  deformation  of  the  nuclear 
equator  and  thus  the  corresponding  geometrical  shadow 
can  be  approximated  by 

e- ft. +!"*«>&*>♦>  .      (3-17) 

This  is  inserted  into  the  expression  for  the 
Praunhofer  amplitude  given  in  (3-15)  which  becomes 

Yc    *** 

•Pre.  =  if  0"*G»<^   e     *  hcM> 
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Equation  (3-18)  is  written  as  a  sum  of  two  integrals, 
one  over  a  disk  of  radius  RQ  and  the  other  over  the 
deformation  distances  as  follows; 

(3-19) 

The  first  of  these  two  integrals  gives  the 
elastic  scattering  component  and  can  be  calculated 

by  noting  the  following  integral  representation  for 

5 

the  Bessel  function. 


■2ro  m 

(3-20) 

Using  (3-20),  it  is  seen  that  the  $  integration  gives 
performing  the  r  integration  gives 


Si.    2.  ZkQo^lt^^ 


(3-22) 


Using  this  result,  the  elastic  angular  distribution 
is  calculated  from  (2-17). 
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Returning  to  the  second  of  the  two  integrals  of 
(3-19),  the  inelastic  scattering  amplitude  will  be 
obtained  where  only  terms  linear  in  the  deformation 
distances  £%tm  will  be  retained.   Integration  of 
the  radial  variable  by  parts  yields. 

where  <A  =  £zk.Sin.M.  Cos<j> 

the  exponential  becomes,  to  the  first  order, 


When  this  is  inserted  into  (3-23),  the  integral 
simplifies  to  the  following. 

J  IM 

o 

To  calculate  the  integral,  the  following  expression 
is  used  , 

S.*n)    c-«).'.'<:^«)!.'    .  .  u   "  , 


itM 


O   foe.  L.-t-M  or>& 


61 


where  (2n)JI  =s  2»4«6  •••(2n) 

and    (2n  +  1)  11   n  1«  3-5  •••(2n  +  1) 
This  is  substituted  in  (3-25),  and  the  order  of 
summation  and  integrations  reversed. 

2xr  (3-27) 

L  +  M  even 
A  consideration  of  the  integral  shows  that  it  can 
be  written  y 

When  m  is  an  integer,  the  Bessel  functions  of.  order 
m  and  -m  are  not  independent  and  are  related  by 
the  expression 

Hence  (3-27)  can  be  written  in  a  more  convenient  form 

*  ^  IMK  JlHl      2- 

The  scattering  amplitude  is  then  multiplied  by 
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v^(  £~  )  and  the  projection  of  this  product  on  the 
final  state  v^(  £;  )  calculated.   The  resulting 
amplitude  becomes 

y    EuM)!(uH)qV'  Tfofagasjwg)  (5"29) 


The  calculation  of  the  matrix  element  requires 

that  a  nuclear  model  be  specified.   In  the  collective 

v 

model,  the  nucleus  is  represented  by  a  deformed 
sphere  of  change.   Excited  states  can  be  characterized 
by  vibrations  and  rotations  of  this  spheroid.   The 
low  lying  states  would  correspond  to  rotations  about 
the  symmety  axis. 

The  collective  model  assumes  the  nucleus  to  be 
in  a  deformed  state  characterized  by  a  multipole 
of  order  L.   The  radius  of  the  nuclear  surface  is 
then 

where  c<  is  the  polar  angle  with  respect  to  the 


symmetry  axis  and  p>-^   is  the  deformation  parameter 

•f*V» 

corresponding  to  the  L   multipole.   By  the  Legenc 
addition  theorem,  this  is  expanded  into  the  space 
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fixed  coordinate  system. 

R-~  Ro[l  +  Pu(^f£x&Y^l  (3-31)' 

where  the  primes  indicate  the  angular  coordinates 
of  the  symmetry  axis.  To  compute  the  value  of  the 
matrix  element,  the  deformation  operators  %f  jjji     mus"k 
be  expressed  in  terms  of  the  variables  used  in  the 
collective  model,   When  (3-31)  is  compared  with 
(3-16),  it  is  seen  that  for  a  particular  L  state, 

If  it  is  first  assumed  that  the  nucleus  is  initially 
in  a  spin- zero  state  and  then  that  the  angular 
dependence  of  the  nuclear  wave  functions  can  be 
expressed  as  spherical  harmonics,  the  calculation 
of  the  matrix. element  is  simply 


<UKl^|00>=A|/: 


The  integral  is  the  orthogonality  integral  for  the 
spherical  harmonics  and  has  the  value  of  unity. 

<S*v^^>-i£y4  (3"55) 
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The  angular  distribution  of  particles 
inelastically  scattered  from  a  nucleus  initially  in 
a  "spin-zero"  state  and  excited  to  a  given  multipole 
L  is  obtained  from  equations  (2-7),  (3-29)  and  (3-33), 


M 

XjfMI(a^slrt§)         '     (3_34) 

This  is  the  basic  result  of  the  Blair-Drozdov  theory. 
For  the  readers  convenience,  the  elastic  and  inelastic 
angular  distributions  for  the  first  four  excited  states 
are  listed  below. 
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The  characteristic  property  of  equation-  (3-34) 
is  that  cross  sections  for  even  values  of  L  are 
calculated  from  a  linear  combination  of  the  squares 
of  even  ordered  Bessel  functions,  while  those  for 
odd  values  of  L  use  Bessel  functions  of  odd  order. 
Beyond  the  first  maximum,  Bessel  functions  rapidly 
approach  their  asymptotic  form, 

t 

T  _^  Sm(3fcE.SH£-ff--fcP>  (3_35) 

Hence  it  is  seen  that  the  angular  distributions  for 
all  even  parity  transitions  have  the  same  shape. 
It  follows,  then,  that  all  odd  parity  transitions 
have  the  same  shape.  The  tv/o,  however,  oscillate 
exactly  180  degrees  out  of  phase  with  each  other. 
It  is  further  noted  that  the  angular  distributions 
for  all  odd  parity  transitions  are  exactly  in  phase 
with  the  elastic  cross  section,  thus  providing  a 
means  for  parity  identification.   This  rule,  known 
as  the  Blair  phase  rule,  is  usually  v/ell  obeyed  in 
experiment. 

Analysis  of  Data 

A  consideration  of  equation  (3-34)  shows  that 
there  are  two  adjustable  parameters  within  the  context 
of  the  theory;  RQ,  a  quantity  of  the  order  of  the 
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nuclear  radius  and/^^,  the  deformation  parameter. 

R0  is  found  by  fitting  equation  (3-22)  to  the  elastic 

cross  section. 

Fits  to  Ar  °(oC,oC)Ar  '  data  are  shov/n  in  figures 
(5-2)  and  (3-3)  for  incident  energies  of  32,8  and 
41.0  Mev(Lab),  respectively,   RQ  was  taken  as  6,45 
Permis  in  both  cases.  For  those  angles  considered, 
the  Blair  theory  predicts  the  oscillatory  nature  of  the 
cross  section,  reproducing  the  maxima  and  minima  very 
well.   It  does,  however,  unlike  the  experimental 
distribution,  have  zeroes.   In  the  back  angles,  the 
magnitude  of  the  theoretical  prediction  drops  off 
more  slowly  than  the  data  and  fits,  in  general,  become 
much  worse. 

It  should  be  noted  that  due  to  the  large 
variations  of  the  elastic  cross  section,  these  figures  . 
show  plots  of  the  ratio  to  Rutherford  scattering, 

Q 

where  the  Rutherford  cross  section  is  given,  by 

With  the  radius  parameter  RQ  determined,  the 
inelastic  distribution  may  be  calculated  to  within 
a  normalization  constant  using  expression  (3-34), 
The  theoretical  distribution  is  then  normalized  to  the 
experimental  cross  section  and y^^is  obtained. 
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There  is  some  concern,  however,  as  to  which  peak 
is  to  be  fitted  since  it  is  expected  that  the  Blair 
formula  will  be  inaccurate  at  small  scattering 
angles  (^20  )  because  of  its  neglect  of  coulomb 
effects  and  at  large  angles  (^L90°)  because  of  various 
small  single  approximations. 

In  the  derivation,  it  was  assumed  that  particles 
having  impact  parameters  less  then  R0  would  be  absorbed, 
and  those  with  larger  impact  parameters  would  pass 
unaffected,  that  is  they  would  not  be  scattered.  For 
charged  projectiles  this  is  an  obvious  oversimplifi- 
cation, as  certainly  two  charged  particles  passing  at 
distances  of  the  order  of  the  nuclear  radius  would 
experience  a  deflection. 

The  apsidal  distance-,  or  distance  of  closest 

approach,  d,  of  a  charged  particle  scattering  from 

9 

the  coulomb  field  of  a  nucleus  is  given  by 


d=f^/^Cscf) 


where  Z  and  /Vare  respectively,  the  charge  numbers 
of  the  target  and  projectile  nuclei,  E  the  center  of 
mass  energy  of  the  incident  projectile  and  ©,  the 
center  of  mass  scattering  angle.  This  is  solved 
for  ©,  giving, 
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©  »  3   Sin"1 


im-v] 


If  the  apsidal  distance  is  taken  to  be  the  order 
of  RQ,  then  for  41  Mev  alpha  particles  incident 

•7LC 

on  Ar   ,  the  scattering  angle  is  of  the  order  of. 
22  degrees.   The  corresponding  result  for  the  32,8  Mev 
alphas  is  27  degrees.   The  "black"  nature  of  the 
nucleus  will  not  allow  impact  parameters  less  then 
RQ  and  thus  the  coulomb  scattering  should  fall  off 
beyond  these  angles.   Hence  normalization  of  the 
cross  sections  to  the  second  peak  seems  justified. 

The  level  scheme   for  the  first  three  low 
lying  states  of  Ar   is  shown  in  figure  (5-4). 
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Pits  to  the  Ar36(cC,0c')Ar   data  are  shown  in 
figures  (3-5)  through  (3-8)  for  32.8  Mev  and  41,0 
Mev  alpha  particles.   It  will  be  noted  that  the 
product  /&^R0   is  more  customarily  quoted  in  the' 
literature,  and  will  be  used  here.  The  first 


71 


excited  state  is  strongly  excited  and  it  seen  to  be 
a  2"*"  state.  Hence  it  is  fitted  using  equation  (3-34) 
with  L  set  equal  to  2,   The  results  for  the  lower 
energy  data  are  shown  in  figure  (3-5)  where  /6gR0  a  .82 
Fermis  and  for  the  higher  energy  data  in  figure  (3-6), 
where  ^gRo  ■■  •7^-  Fermis.   It  is  interesting  to  note 
how  well  the  Blair  phase  rule  is  obeyed,   A  comparison 
of  the  L  =*  2  data  with  the  corresponding  elastic  data 
shows  them  to  oscillate  exactly  180  degrees  out  of 
phase. 

The  second  excited  state  is  strongly  excited, 
with  its  spin  known  to  be  either  2  or  3,  Upon 
comparison  of  the  experimental  distributions  shown 
in  figures  (3-7)  and  (3-8),  with  the  corresponding 
elastic  data,  they  are  seen  to  oscillate  in  phase. 
If  the  Blair  phase  rule  is  correct,  then  the  second 
excited  state  can  only  be  a  3"  state  and  thus  its  spin 
and  parity  are  determined.   The  angular  distributions 
are  calculated  using  (3-32)  with  L  set  equal  to  3. 
These  results  are  also  shown  in  figures  (3-7)  and 
(3-8)  where  ^Ro  s   •67  P©™!-3  ?or   the  lower  energy 
data,  and  .56  Fermis  for  the  higher  energy  data. 

The  general  appearance  of  the  inelastic  angular 
distributions  is  similar  to  their  elastic  counterparts. 
The  oscillatory  nature  is  predicted  and  the  position 
of  the  maxima  and  minima  are  reproduced  well. 
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As  in  the  elastic  case,  the  inelastic  distribution 
predicts  zero  minima  and  fares  poorly  in  the  back 
angles. 

The  Distorted  Wave  Born  Approximation 

,  This  and  the  remaining  text  follow  the  treatments 

11  12 

given  by  Austern   and  McCarthy, 

If  it  is  assumed  that  the  average  interaction 

potential  K(r)  of  the  colliding  particle  is  primarily 

responsible  for  determining  the  optical  wave  functions 

X.(k,r),  and  if  k*  tfc  \cr ,  then  the  potential  V"(r,  »  ) 

may  be  treated  as  a  perturbation.   Consider  the 

matrix  element  T^i  giving  the  inelastic  component  of 

the  scattering, 

T«»  <^XW|V)^>  (3-36) 

The  wave  function  ^  is  then  approximated  by  the 
outgoing  optical  wave  function  that  was  used  in 
(2-41)  to  compute  the  elastic  scattering  component. 


If  *  V^X^i:)  .        ,   (3-37) 


where  the  approximation  of   (3-1)   has  been  applied. 
With  this  expression  f or  ^£> ,   equation   (3-36)   becomes 
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TF1    -Jfrjrt&lYltyXS  (3-38) 

The  width  factor  ^vf|Vf lvA  Influences  the  magnitude 
of  the  cross  section  and  depends  on  a  detailed 
knowledge  of  both  the  reaction  processes  and  the 
nuclear  wave  functions.   The  distorted  wave  functions 
/C  (k,r),  are  primarily  responsible  for  determining 
the  angular  distribution. 

The  potential  V^Jr,!?),  assumed  a  scalar  function 
of  position,  Is  expanded  In  multipoles  of  the  vector  r. 
This  expansion  becomes 

where  it  is  necessary  for  the  ensuing  calculations 
that  VTr.  behave  like  the  spherical  harmonics  under 
rotations  of  the  coordinate  system  and  have  parity 
(-1).  The  factor  ^vi.|Vf|vi\  can  then  be  written. 

<WIS&  "ZKW^VJY*  (3-40) 

l-M 

The  wave  functions  v^  and  v^  are  assumed  to 
have  their  angular  dependence  described  by  the 
spherical  harmonics*  They  will,  in  general,  be  a 
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function  of  their  angular  momentum  Ji\   and  its 
projection  ISIS.   A  fundamental  property  of  the  matrix 
elements  of  these  spherical  harmonics  and  functions 
v/hich  transform  in  a  like  manner  under  a  coordinate 
rotation,  is  the  Wigner-Echart  theorem.  ^ 

Applied  to  the  above  "width"  factor,  it  is 

<vWviHJW=  QcLMxMU,nXVlvM>   (3"41) 

J.  and  J^  are  the  initial  and  final  nuclear  spins 
with  Mjt  and  AU  their  respective  projections.   The 
factor  <vj  }]  VL  \l   Vj A  is  the  "reduced  width"  and 
(J*L  M±M  |  JriMf)  a  Clebsch-Gordan  coefficient.   The 
reduced  width  is  just  the  width  factor  expressed  in 
a  diagonalized  representation.   Using  (3-41),  equation 
(3-40)  becomes 

<v„ivivi>  -Z"C4L»4MiJW«gi\(i*>x£  (3~42) 

uii 

■f*V» 

The  interpretation  being  that  the  L   multipole  in 
equation  (3-39)  transfers  an  amount  of  angular  momentum 
L  to  the  nucleus.   The  Clebsch-Gordan  coefficients 
are  so  designed  as  to  limit  L  to  the  range 

Jix-Jfl  i  L£|J_r+JPJ 

and  thus  insure  conservation  of  angular  momentum. 
The  "reduced  width"  factor  is  a  function  of  only  the 


/ 
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r  variable  and  It  becomes  convenient  to  write  it 
as  the  product  c^'  a  "strength  factor"  and  a  "form 
-actor"  which  gi"  3  its  radial  dependence. 

<Y>P!!VJ!V ;  ■  ■-  ALFL(.f)  0-43) 

The  strength  factor  A^  would  give  the  magnitude  of 
the  reaction  and  requires  a  specific  knowledge  of  both 
the  nuclear  wave  functions  and  the  nature  of  the 
reaction.   The  form  factor,  however,  need  not  be 
known  exactly,  and  meaningful  results  are  obtained 
from  knowing  its  general  trend  as  will  be  seen  later. 

Upon  substitution  of  expressions  (3-42)  and 
(5-43)  into  (3-38)  the  matrix  element  becomes 


V  - 


where  (3-44) 

The  cross  section  is  then  obtained  from  expression 
(2-110), 

AV   UM  (3-45) 

It  will  be  recalled  that  the  summation  sign  referred 
to  a  sum  over  UL>   and  an  average  of  M^.  For  an 
initial  state  having  spin  Ji#  there  are  (2J.4 J)  values 


80 
of  M^,   hence    (3-45)   becomes 


tfc        Mp     UK     fcV  (3-46) 

Noting  the  following  Clebsch-Gordan  symmetry 

14 
relationship, 

the  cross  section  can  be  written 


u  .iv^ry^ /teg)  frjr+tf   V"  S"i  a  »*i<v  l* 


The  Clebsch-Gordan  coefficients  exhibit  the  following 
closure  relationship • 


Hence  the  cross  section  expression  simplifies  and 
is  written 

~„I!£2^/'!£iA/^,lf^^7_'i  7T  I2   (3-50) 

If  the  spin  of  the  initial  state  is  zero,  then  the 
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final  state  spin  is  just  the  angular  momentum 
transferred  in  the  reaction.   It  should  be  noticed 
that  different  values  of  the  momentum  transfer  L 
do  not  interfere  in  determining  the  reaction. 

The  problem  of  evaluating  the  cross  section 
has  now  reduced  to  that  of  calculating  the  matrix 
element  £.  L--.   Her  i       Lt  .uiias  of  the  optical  wave 
functions  are  required.   These  may  be  approximations   ' 
or  with  complete  rigor,  the  spherical  harmonic 
expansions  derived  in  chapter  II  may  be  used.  The 
following  sections  discuss  the  nature  of  possible 
approximations • 

Plane  Waves 

While  it  may  seem  paradoxical  to  include  a 
section  with  this  title  in  a  discussion  of  the 
distorted  wave  Born  approximation,  the  use  of  plane 
waves  is  instructive  as  it  gives  insight  into  the 
nature  of  the  distortion  necessary  to  reproduce  the 
form  of  the  experimental  angular  distributions. 

It  is  convenient  to  make  the  identification 

with  the  matrix  element  Cl^I  given  in  (3-45) 
then  becoming 
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The  quantity  /?(r)  eiS^r'  is  a  probability  density 

and  may  be  thought  of  as  describing  the  properties 

of  a  probe  which  measures  the  evolution  of  the  nuclear 

state  during  the  process  of  inelastic  scattering. 

To  clarify  this  idea,  consider  the  case  where  the 

optical  wave  functions  have  been  replaced  by  plane 

waves. 

where  K  =r  kj.  -  kf  and  pQ   is  a  constant.  Here  the 
probe  measures  the  difference  in  the  wave  number 
vectors  for  the  initial  and  final  states  and  the 
matrix  element  2TlJ£  is  proportional  to  the  probability 
that  the  initial  and  final  states  differ  by  momentum 
KK.   Since  the  inelastic  scattering  cross  section  is 
obtained  from  this  matrix  element,  the  use  of  plane 
waves  implies  that  the  entire  momentum  transferred 
goes  into  the  inelastic  component  of  the  scattering. 

If  /?(rj  is  anisotropic  in  space,  the  probe  tries 
to  make  a  simultaneous  measurement  of  its  position 
and  the  momentum  transfer.   By  the  uncertainty  principle, 
the  resolution  of  the  momentum  transfer  measurement 
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18 
must  then  be  reduced.   In  But  lei*1  s   theory  of  nuclear 

reactions,  it  is  assumed  that  the  reaction  is  confined 

to  the  surface.   It  would  then  seem  that  the  probe 

could  obtain  no  information  whatsoever  about  the 

radial  component  of  the  momentum  transfer.   On  the 

other  hand,  since  no  angular  localization  is  assumed 

on  the  nuclear  surface,  the  resolution  of  the  angular 

momentum  transfer  should  be  a  maximum.   This  is  seen 

in  the  following  calculation. 

Setting  the  constant  /Q 0  equal  to  unity,  the 

matrix  element  TSjjj  D©c°^©s 

where  P(r)  becomes  a  Dirac  delta  function  limiting 
the  reaction  to  a  spherical  shell  of  radius  R0. 
Using  equation  (2-18),  the  plane  wave  is  expanded 
in  a  series  of  spherical  harmonics  and  (3-54)  becomes 

The  integral  in  (3-55)  is  the  orthogonality  integral 
for  the  spherical  harmonics  where 

JYmYwcJSI   -  $*!&-*'  (3-56) 


J 
J 
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The  Kronecker  delta  functions  are  zero  unless  L1  and 
M1  equal  L  and  M  respectively,  hence  only  the  LM 
term  remains. 
Equation  (3-55)  becomes 

where  the  primes  indicate  the  angular  coordinates 
of  the  symmetry  axis.   As  it  will  be  assumed  that  the 
nucleus  is  initially  in  a  spin  zero  state,  the  cross 
section  for  a  given  final  state  L  is  given  by  (3-50) 
to  be  proportional  to 

r~ 

oi *  \iMWxtifr*)Y&  (3"58) 

The  Legendre  addition  theorem  is  used  to  convert  the 
sum  of  spherical  harmonics  into  a  Legendre  polynomial 
of  argument  unity  which  has  magnitude  unity  for  all 
values  of  L.   Equation  (3-58)  becomes 

a:*  ji(^)  l3"59) 

Beyond  their  first  maximum,  these  spherical 
Bessel  functions  rapidly  approach  their  asymptotic 
form  given  in  equation  (2-19),  and  thus  in  this  region, 
the  cross  section  can  be  written, 


v 
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01 U^^  (3"50) 


Hence  this  behaves  like  the  Blair  formula.   It 
appears  then,  that  the  sharpness  of  the  maxima  and  the 
fact  that  the  minima  actually  reach  zero  can  be 
associated  with  the  maximum  resolution  of  the  angular 
momentum  transfer.   If  this  interpretation  is  correct, 
it  should  be  possible  to  spoil  the  resolution  by 
causing  /?(r)  to  be  non- uniform  over  the  nuclear 
surface. 

We  proceed  by  localizing  the  reaction  to  part 
of  the  nuclear  surface.   As  alpha  particles  are 
strongly  absorbed,  those  incident  on  the  nuclear 
surface  will  have  their  wave  funcions  considerably 
reduced  on  the  shadow  side  of  the  nucleus  as  compared 
v/ith  the  illuminated  side.  The  same  effect  would  be 
true  of  the  scattered  particles.   The  overlap  of  their 
wave  functions  would  be  greatest  on  the  side  of  the 
nucleus  opposite  to  the  direction  of  momentum  transfer 
as  is  shown  in  figure  (3-9).   It  is  assumed  that 
that  part  of  the  surface  axially  symmetric  about  the 
direction  of  momentum  transfer  and  cut  off  by  a  cone 
of  semiangle  o( contributes  uniformly  to  the  matrix 
element,  and  the  remaining  part  of  the  surface 
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FIGURE  (3-9) 

contributes  nothing  at  all.   The  matrix  element  2*tm 

is  calculated  using  the  expression  in  (3-54)  where 

the  plane  wave  has  been  expanded  in  a  series  of 

Legendre  polynomials.   After  the  radial  integration, 

(3-54)  becomes 

/ 


(3-61) 


where  the  axis  of  symmetry  is  taken  along  the  direction 
of  momentum  transfer.  The  Legendre  polynomials  and 
spherical  harmonics  are  related  through  the  expression 


#*•>■  j!  x#*> 


(3-62) 


and  as  the  azimuthal  integration  has  been  unaffected, 
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it  is  seen  by  (3-56)  that  the  integral  is  non-zero 
only  when  M  equals  zero.   Using  this  result,  the 
matrix  element  Ltm   becomes 


A'm,LA% 


(3-63) 


The  integrand  is  expanded  in  a  series  of 
Legendre  polynomials. 


SL  (3-64) 

- p*  (g,  Xf*> 

The  expansion  coefricient  is  obtained  by  multiplying 
both  sides  of  (3-64)  by  //  (Oj(j))  d&and  integrating 
over  4TC  radians.   The  orthogonality  of  the  spherical 
harmonics  determines  the  coefficient  Cj^1  j£   on  the 
right  thus  giving  the  expression 

<£&'  *Jf %&%><**•  '  (3"65) 


The  following  relation  involving  Clebsch- (Jordan 

19 
coefficients  is  noted, 

J  (3-66) 
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Using  this  relation,  equation  (3-55)  becomes 

CuV*    £&+,»„  J*  <LLooUo)  (3"67) 

The  expansion  in  (3-64)  is  then  obtained. 

^  ^  m  teg^^  (3-68) 

JL-.o 
The  matrix  element  *£j„  becomes 


l'ji        *  g 
X J  %  &>s&)  Strode 


(3-69) 


Hence  for  a  given  angular  momentum  transfer  L,  the 
matrix  element  is  proportional  to  a  sum  of  spherical 
Bessel  functions  of  different  orders.   In  general 
this  will  be  complex  and  it  is  therefore  highly 
improbable  that  2Ttvj  will  go  to  zero  for  any  real 
value  of  K.   The  angular  resolution  associated  with 
zero  minima  has  been  spoiled.   If  &  =  TC,   the  entire 
surface  then  contributes  uniformly.   The  integral 
in  (3-69)  becomes  the  orthogonality  integral  for  the 
Legendre  polynomials,  giving  zero  unless  <KS0* 
With  A  «  0,  the  Clebsch- (Jordan  coefficient  will  be 
non-zero  only  when  L1  s  L,  hence  it  is  seen  that  the 
result  reduces  to  that  in  equation  (3-59). 
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An  interesting  interpretation  of  (3-69)  is 
that  for  total  angular  momentum  transfer  L,  L1 
represents  that  going  into  the  inelastic  component 
of  the  scattering  as  it  results  from  the  L1   partial 
v/ave  in  the  expansion  of  a  plane  wave,  and  k.   that  which 
has  gone  into  other  components,  such  as  elastic 
scattering  or  absorption  from  the  beam. 

Approximate  Optical  Wave  Functions 

20  21 
Based  on  a  study  of  optical  model  wave  functions,   » 

22 

reasonable  approximations   for  the  incoming  and 

outgoing  wave  functions  are  found  to  be 

2    A 


%£  *  Mr)  0    ~" 

where  Jfis  referred  to  as  the  anisotropy  parameter 
and  is  a  measure  of  the  non  uniformity  of  the  nuclear 
surface.  When  these  expressions  are  used  in  equations 
(3-51)  and  (3-52),  the  matrix  element  *£"m  is  seen 
to  be 

Where  as  usual,  jC  ■»  k^  -  k£ .   Let, 
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if-**  a  <£-&,) 


t  -  <£•£) 


'Ac. 


(3-72) 


and  expand  the  resulting  plane  wave  in  a  series  of 
spherical  harmonics.  The  angle  integrations  were 
computed  in  equations  (3-55)  and  (3-56)  and  they 
result  in  the  following  expression  for  Ltm* 


Xh  -  ^il^^4m^M^AM./^    (3-73) 

a 

The  cross  section  is  calculated  using  expression 
(3-50),  where  the  spin  of  the  initial  state  is  assumed 
to  be  zero. 


OL*Z)£fo#Ytiv& 


-/.■ 


Using  the  Legendre  summation  formula, 

mm  mm 

where  P  is  the  angle  between  ^  andj£  . 
CosCP)can  be  written 


2.  (3-74) 


(3-75) 


c*s<r)  - 


(3-76) 
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IT  the   adiabatic   approximation  ki  ft  kf   «  k,    is 
employed,^  can  be  written 


£**   Z(K  +  cJt)S/n@     &  (3-77) 


where  &  is  the  angle  between  k*    and  k^.  and  K,  is  a 
unit  vector  along  the  direction  of  momentum  transfer, 
(3-75)  is  inserted  in  (3-74)  and  it  can  be  shown  that 
Cos(  F    )   =2.  1.   Hence  PL(Cos  P  )   ^  1  and  the  cross 
section  for  a  given  value  of  L  is  proportional  to 


01  oc 


z 


(3-78) 


a 

It  now  remains  to  make  some  approximation  for 
the  radial  integral.   Since  alpha  particles  are 
strongly  absorbed  by  the  nucleus,  I  A(r)  |  ''  will 
fall  off  for  r  less  than  the  nuclear  radius.   On  the 
other  hand,  P(r)  being  associated  with  the  nuclear 
wave  function,  will  fall  off  for  r  greater  than  the 
nuclear  radius.   Assuming  that  P(r)  is  slowly  varying 
over  the  mean  path  of  the  alpha  particle,  the  radial 
integration  should  essentially  depend  on  the  radius  RQ 
and  the  thickness  of  the  surface.   A  simple  empirical 
choice  which  involves  only  these  two  parameters  is 
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The  final  form  for  the  differential  cross  section 
is  therefore 


<t* 


.'£►-' 


(3-80) 


6 

With  the  adiabatic  approximation  for  gC  given  in 
(3-77),  (3-80)  can  be  written  in  a  simpler  form. 


dec 


7    />•-£•  \* 

f 


The  corresponding  surface  interaction  result  is 


<TL  cC 


jut2tkRo+CY)SbJl 


(3-82) 


In  the  discussion  resulting  in  equation  (3-69)  it 
was  found  that  by  introducing  an  angular  ani sot rophy 
into  the  reaction,  the  resolution  for  a  measurement 
of  the  angular  momentum  transferred  was  at  least 
partially  spoiled. 

It  was  seen  that  the  matrix  element,  T^i^,   was 
no  longer  proportional  to  a  single  spherical  Bess el 
function  of  order  L,  but  consisted  of  a  sum  of  Ajj(t££) 
The  result 'in  (3-82)  is  but  a  single  spherical  Bessel 
function.  However  its  argument  is  complex.  That  this 
actually  corresponds  to  (3-69)  is  shown  by  expanding 
(3-82)  in  a  series  of  jLf  (2KR0Sin^. 
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Consider  the  plane  wave  expansion,  where  the 
magnitude  of  K  RQ  a  X  **y 

Z(jt4>y)CosG 

This   can  also  be  written 


xfa (*)  fa  Ly)  Pu,  Ceo <£>)  gccose) 


(3-84) 


The  two  series  are  equated  and  both  sides  multiplied 
by  PL(Cos9)  Sin9d9.   When  the  result  is  integrated 
from  o  to  TC  ,   orthogonality  eliminates  all  but  the  - 
L'  '  term  in  (3-83)  giving  the  following  expression 

-itOc+y)  -  Y^A'L  fasten)  <U' fa  AA  (y) 

U/  (3-85) 

t   /       , ./A  f 

The  integral  is  written  in  terms  of  spherical  harmonics 
and  both  sides  are  multiplied  by  d<j)  and  integrated 
from  0  to  ZTC 

The  integral  is  given  in  (3-66),  and  with  the  proper 
substitutions  (3-86)  becomes 
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jutx+y)  =2  tl**~l &*'+')("'*<> i£*)xpto<u(y)  (3-8?) 


Letting, 


It  is  seen  that- 


3-88) 


which  is  the  desired  result.  The  parameter  Jr 
evidently  has  the  same  effect  on  reducing  the  resolution 
of  an  angular  momentum  transfer  measurement  as  the 
introduction  of  an  angular  anisotropy  into  the  reaction. 

Results 

Optical  model  fits  to  the  Ar36  (  o< ,  c<)  Ar3  data 
are  shown,  for  the  lower  energy  case  in  figure  (3-10), 
and  in  (3-11)  for  the  higher  energy  case.  These  are 
included  as  the  approximations  in  (3-68)  are  based  on 
a  study  of  the  optical  wave  functions . and  thus  the 
ability  of  the  optical  model  to  predict  the  elastic 
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distributions  should  indicate  the  validity  of  these 
approximations  for  the  data  considered.  The  maxima 
and  minima  are  reproduced  very  well,  both  in  position 
as  in  the  Blair  formula,  and  in  magnitude.   Although 
not  indicated,  the  optical  model  does  well  in 
predicting  the  behavior  of  the  back  angles,   A  radius 
of  5,6  Permis  was  used  in  both  cases. 

Pits  to  the  Ar36((X,  OU)Ar36^ata  obtained  using 
(3-81)  and  the  appropriate  choice  of  L  are  shown  in 
figures  (3-12)  through  (3-15),  A  single  set  of  the 
parameters,  A   and  V",  were  used  to  fit  these  data. 
They  were  A*  0*88  and  V—  0,90  as  was  suggested  by 
McCarthy  and  Pursey  in  the  text  of  their  article „ . 
Similarly  a  single  value  of  the  radius  parameter  R^ 
was  used.   This  was/ found  to  be  7,0  Permis,   It  is 
seen  that  the  angular  predictions  of.  the  maxima  and 
minima  are  comparable  to  the  Blair  theory.  There  is, 
however,  improved  agreement  in  the  relative  magnitudes 
of  the  theoretical  and  experimental  distributions.. 
The  non-zero  minima  associated  with  the  poorer 
resolution  of  the  angular  momentum  transfer  measurement 
is  clearly  shown.  A  systematic  study  of  the  parameters 
/  and  Y  was  not  made.  The  effect  of  Y»   however, 
can  be  inferred  from  the  previous  study.  Setting  Y 
equal  to  zero  removes  the  anisotropy  and  thus  the 
minima  will  again  reach  zero.  The  interpretation 


98 

of  the  series  in  (3-88)  was  that  the  entire  angular 
momentum  transfer  Lh  no  longer  went  into  the  inelastic 
component  of  the  scattering  but  only  an  amount  L*fi« 
The  rest  was  absorbed  out  of  the  beam  or  went  into 
the  elastic  scattering  component.  Hence  it  is  expected 
that  Y"  would  also  affect  the  rate  at  which  the  cross 
section  decreased  with  angle* 
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APPENDIX  A 

THE  CONFLUENT  HYPERGEOMETRIC  EQUATION 

Hummer1 s  equation  is  the  hype rgeome trie  equation 
often  encountered  by  physicists  in  the  solution  of 
Schrodinger • s  equation. 
It  has  the  form 

(A-l) 


2d4  +  (b-z)d!±-aS=0 


Solutions  of  this  equation  are  known  as  confluent 
hypergeometric  functions.   This  appendix  will  present 
only  enough  of  the  theory  to  provide  an  understanding 
of  chapter  II.  The  reader  is  referred  to  references 
£~\J   an<^  C^J  £°r  a.  more  complete  development. 

The  Series  Solution 

It  is  assumed  that  equation  (A-l)  has  a  solution 
of  the  form 

OO 


M  =  Zc»z 


Ic+n 


This  is  inserted  into  the  differential  equation 
resulting  in  a  polynomial  in  z.  For  y  to  be  a  solution 
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for  all  values  of  z,  the  coefficients  of  all  powers 
of  z  must  be  identically  zero.   With  n  set  equal  to 
zero,  the  coefficient  of  the  lowest  remaining  power 
of  z  yields  the  indicial  equation. 
Solving  this,  gives 

k  =  0,  1-b  (A-2) 

The  first  of  these  gives  the  desired  solution  which  is 
regular  at  the  origin.   The  following  recurrence 
relationship  for  the  coefficients  is  then  obtained, 
assuming  that  c0  is  non-zero. 

CMl  =    &±g       cn  u-3) 

With  this,  all  the  coefficients  can  be  generated 
to  within  an  arbitrary  constant  c  and  setting  this 
equal  to  unity  yields  the  following  series 

A  Contour  Integral  Representation 
The  expression  (A-4)  is  re-written 

Ff*iu*)-       r(b)      yr^rdr-tiz" 


rc«&rar 
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where  B(u,v)  is  the  beta  function  and  is  defined 


r(u+v) 


*->. 


tu-'o-t)'dt- 


Jo 


A  contour  integral  representation  of  the  beta  function 
can  be  shown  to  be 


bm=  o-^ydturdt 


(A-6) 


where  the  contour  C-  is  described  in  figure  (A-l)  and 
the  phase  defined,  so  that  along  the  segment  from 
1  to  2, 


arg 

[t)   m  arg   (1-t)  =  0 
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r 
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1 

FIGURE  (A-l) 

For  the  expression  to  take  on  the  same  value  after  a 
complete  circuit  of  the  contour,  it  is  seen  that  u  +  v 
must  equal  an  integer.  Equation  (A-6)  is  inserted 
into  (A- 5)  giving 
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Interchanging  the  order  of  summation  and  integration, 
where  it  is  seen  that  b  must  be  an  integer  to  satisfy 
the  requirement  that  (A- 6)  be  single  valued  and  thus 
(A-7)  be  a  solution  to  (A-l),  this  becomes 


'  rMr(b-o.)      j     rt' 

(A-7) 

rc^rCb-a)  J 

S3 

When  the  integrand  of  (A-7)  is  substituted  back 

into  the  differential  equation,  the  result  can  be 
expressed  as  a  perfect  differential  and  thus  any  contour 
which  has  the  same  value  at  its  end  points  for  v/hich 
the  integral  converges  leads  to  a  solution.   In 
particular,  two  of  these  will  be  considered  yielding 
irregular  solutions  of  (A-l) • 

Let      ,IT(a/b|2)^V;^/b|a)4^aibi2)      u.8) 

where 

MAlfchQ.  r^^g^rf/H<)^-    (a-9) 

*C/9CZ 

The  contours  C^  and  C2  are  described  in  figure 

(A-2) 


FIGURE  (A- 2) 
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where  for  convergence  of  the  integrals,  ■ 

is  the  argument  of  the  end  points  at  infinity. 
Consider  the  expression  (A-9) .  By  convention, 
the  substitution  zt  =  -u  is  made  and 


W/alw*) 


vm-?^L  ruitUta'k 


r<b)n  (*>-«•)  J 


(A- 10) 


where  the  contour  in  the  U-plane  is  shown  in  figure 
(A-3) 


Crt       wyCtc)-* 


FIGURE  (A-3) 
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Here,  arg  (u)  ss  0  along  the  bottom  portion  of  the  curve 
along  the  real  axis  as  was  defined  in  figure  (A-l) . 
With  this  contour, 

Letting  first  t  =  t  «fr  1  and  then  zt  as  -u,  an  expression 
for  W2(alb|z)  is  obtained 


The  contour  C&  is  shown  in  figure  (A-4) 

A 


(A- 11) 
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FIGURE  (A-4) 

where  arg(u)  a;0  on  the  top  portion  of  the  curve  along 
the  real  axis,  again  to  be  consistent  with  the  phase 
convention  in  figure  (A-l),   To  obtain  convergence 

-7%  <  a.*y  2  <  K/2. 

Asymptotic  Expansion 

Asymptotic  series  expressions  will  be  obtained 
for  the  irregular  solutions  W\  and  W2  and  from  these 


Ill 


and  (A-8) ,  an  asymptotic  expression  Tor  ^P^Catblz) 
v/ill  be  found. 

Consider  expression  (A-10)  and  its  associated 
contour.   The  binomial  expansion  is  used  to  obtain 


h*o         J 


This  is  a  series  of  negative  powers  of  z  and  is  thus 
an  expansion  for  large  z.   This  is  inserted  into 
(A-10)  and  the  order  of  summation  and  integration 
reversed. 


£«~ 


1: 


The  remaining  integral  and  its  associated  contour 

C4  given  in  figure  (A-3)  can  be  shown  to  be  a 

contour  integral  representation  of  the  gamma  function. 

Equation  (A-ll)  becomes 

This  is  put  in  a  form  where  (-z)  is  the  variable 
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of  the  series  with  the  following  identity  which  can  be 
verified  by  induction. 

_£t±L  =  ^.r(o.^0  u.13) 

where  n  is  an  integer.   Equation  (A- 12)  becomes 

In  expression  (A-ll),  the  binomial  theorem  is  used 
to  obtain 

A  \i'K  f*  PW  ^ 

Equation  (A-ll)  becomes 

As  it  was  required  that  b  be  an  integer,  this  can 
be  written 


(A-15) 
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p-<o 


v-Al 


The  remaining  integral  and  its  related  contour  can  be 
shown  to  be  a  contour  integral  representation  of  the 
gamma  function 

Equation  (A- 16)  becomes 

w>i^j  „  cfl^+j:  caicfefial  ^""  (A-i7) 

Letting  a  -^- a  in  (A- 13),  equation  (A- 17)  can  be 
transformed  in  a  manner  similar  to  (A- 14) 


VLfelbtf  *.  -S&  ^^T^-i-qjK^M^  u.18) 

For  very  large  z,  only  the  first  term  of  the 
expansion  is  retained,  thus  equation  (A-8)  gives  for 
the  regular  asymptotic  solution 
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>fT(a/b|3)  *  ^JSL(^%  ^   gg-*°  (A-19) 


where 


-Zffc  <  Q-^^  <  ^/^ 


In  the  development  used  in  chapter  II,  z  takes  on 
pure  imaginary  values,  '  > ♦ 

arg(z)  =  ±  *?/& 
In  order  to  use  expression  (A-19)  a  limiting  procedure 
will  have  to  be  conceived  where  arg(z)  is  allowed  to 
approach  f/2  or  -^2  *n  *ke  direction  of  positively 
or  negatively  increasing  angles,  respectively.   It  will 
be  assumed  that  this  is  done  in  all  expressions  therein, 
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